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Abstract 

The credit crisis and the ongoing European sovereign debt crisis have highhghted the 
native form of credit risk, namely the counterparty risk. The related Credit Valuation 
Adjustment (CVA), Debt Valuation Adjustment (DVA), Liquidity Valuation Adjust- 
ment (LVA) and Replacement Cost (RC) issues, jointly referred to in this paper as 
Total Valuation Adjustment (TVA), have been thoroughly investigated in the theoreti- 
cal papers |Crl2a| and |Crl2b| . The present work provides an executive summary and 
numerical companion to these papers, through which the TVA pricing problem can be 
reduced to Markovian pre-default TVA BSDEs. The first step consists in the counter- 
party clean valuation of a portfolio of contracts, which is the valuation in a hypothetical 
situation where the two parties would be risk-free and funded at a risk-free rate. In 
the second step, the TVA is obtained as the value of an option on the counterparty 
clean value process called Contingent Credit Default Swap (CCDS). Numerical results 
are presented for interest rate swaps in the Vasicek, as well as in the inverse Gaussian 
Hull- White short rate model, also allowing one to assess the related model risk issue. 

Keywords: Counterparty risk. Credit valuation adjustment (CVA), Debt valuation ad- 
justment (DVA), Liquidity valuation adjustment (LVA), Replacement cost (RC), Backward 
stochastic differential equation (BSDE), Levy process. Interest rate swap. 

1 Introduction 

The credit crisis and the ongoing European sovereign debt crisis have highlighted the native 
form of credit risk, namely the counterparty risk. This is the risk of non-payment of promised 
cash-flows due to the default of the counterparty in a bilateral OTC derivative transaction. 
The basic counterparty risk mitigation tool is a Credit Support Annex (CSA) specifying a 
valuation scheme of a portfolio of contracts at the default time of a party. In particular, the 
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netting rules which will be applied to the portfolio are specified, as well as a collateralization 
scheme similar to margin accounts in futures contracts. 

By extension counterparty risk is also the volatility of the price of this risk, this price 
being known as Credit Valuation Adjustment (CVA), see |PR05| . |BMP13| . [BBCH13]. 
Moreover, as banks themselves have become risky, counterparty risk must be understood in 
a bilateral perspective, where the counterparty risk of the two parties are jointly accounted 
for in the modeling. Thus, in addition to the CVA, a Debt Valuation Adjustment (DVA) 
must be considered. In this context the classical assumption of a locally risk- free asset 
which is used for financing purposes of the bank is not sustainable anymore. This raises 
the companion issue of proper nonlinear accounting of the funding costs of a position, and 
a corresponding Liquidity Valuation Adjustment (LVA). The last related issue is that of 
Replacement Cost (RC) corresponding to the fact that at the default time of a party the 
contract is valued by the liquidator according to a CSA valuation scheme which can fail to 
reflect the actual value of the contract at that time. 

The CVA/DVA/LVA/RC intricacy issues, jointly referred to henceforth as TVA (for 
Total Valuation Adjustment), were thoroughly investigated in the theoretical papers |Crl2a| . 
|Crl2b| . [BCPPlT] and |BK11] . The present work is a numerical companion paper to 
|Crl2aj and |Crl2bj . Section [2] provid es an executive summary of the two papers. Section [s] 
describes various CSA specifications. In Sections |4] and [s] we present clean valuation (clean 
of counterparty risk and funding costs) and TVA computations in two simple models for 
interest rate derivatives. We show in a series of practical examples how CVA, DVA, LVA 
and RC can be computed in various situations, and we also assess the related model risk 
issue. 



2 TVA Representations 
2.1 Setup 

We consider a netted portfolio of OTC derivatives between two defaultable counterparties, 
generically referred to as the "contract between the bank and her counterparty". The 
counterparty is most commonly a bank as well. Counterparty risk and funding cash-flows 
can only be considered at a netted and global level, outside the scope of different business 
desks. Consequently, the price 11 of the contract must be computed as a difference between 
the clean price P provided by the relevant business desk and a correction Q computed 
by the central TVA desk. By price we mean here the cost for the bank of margining, 
hedging and funding ("cost of hedging" for short). By clean price we mean the price of 
the contract computed without taking into account counterparty risk and excess funding 
costs. Symmetrical considerations apply to the counterparty, but with non-symmetrical 
data in the sense of hedging positions and funding conditions. As a consequence and due 
to nonlinearities in the funding costs, the prices (costs of hedging) are not the same for the 
two parties. For clarity we focus on the bank's price in the sequel. 

We denote by T the time horizon of the contract with promised dividends dDt from 
the bank to her counterparty. Both parties are defaultable, with respective default times 
denoted by 9 and 9. This results in an effective dividend stream dCt = JtdDt, where 
Jt = lt<r with r = 9 A 9. One denotes by f = r A T the effective time horizon of 
the contract as there are no cash-flows after f. The case of unilateral counterparty risk 
(from the perspective of the bank) can be recovered by letting 9 = oo. After having sold 
the contract to the investor at time 0, the bank sets-up a collateralization, hedging and 
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funding portfolio ("hedging portfolio" for short). We call an external funder of the bank 
(or funder for short) a generic third-party, possibly composed in practice of several entities 
or devices, insuring funding of the position of the bank. This funder, assumed default-free 
for simplicity, thus plays the role of "lender /borrower of last resort" after exhaustion of the 
internal sources of funding provided to the bank via the dividend and funding gains on her 
hedge, or via the remuneration of the margin amount. 

The full model filtration is given as G = F VM^ VM^ where F is a reference filtration 
and n1 = a{9 A t), Tif = a(e A t). A probability space (fi,^T,lP), where P is some risk- 
neutral pricing measure, is fixed throughout. The meaning of a risk-neutral pricing measure 
in this context, with different funding rates in particular (see |Crl2aj ). will be specified by 
martingale conditions introduced below in the form of suitable pricing backward stochastic 
differential equations (BSDEs); see [EPQ97| for the seminal reference in finance. 

Remark 2.1 Even though it will not appear explicitly in this paper, a pricing measure 
must also be such that the gain processes on the hedging assets follow martingales, see 
|(^rl2aj . 

Moreover, we assume that F is immersed into G in the sense that an F-martingale stopped 
at r is a G-martingale. 

Remark 2.2 As discussed in |Crl2b| . this basic assumption precludes major wrong- way 
risk effects such as the ones which occur with counterparty risk on credit derivatives. In 
particular, under these assumptions an F-adapted cadlag process cannot jump at r. We 
refer to |Crl2cj for an extension of the general methodology of this paper beyond the above 
immersion setup to incorporate wrong- way risk when necessary. 

2.2 Data 

We denote by rt an OIS rate, where OIS stands for an Overnight Indexed Swap, the best 
market proxy of a risk-free rate. By ft = + jt we denote the credit-risk adjusted rate, 
where jt is the F-hazard intensity of r, which is assumed to exist. Let (3t = exp(— rsds) 
and = exp(— rsds) stand for the corresponding discount factors. Furthermore, lEj and 
Et stand for the conditional expectations given Qt and J-j, respectively. The clean value 
process Pt of the contract with promised dividends dDt is defined as, for t £ [0,f], 

PtPt ■= it MDs^ ^^'^[ + 

by immersion of F into G. 

Effectively, promised dividends dDt stop being paid at r (if r < T), at which time 
the last terminal cash-flow R paid by the bank closes out her position. Let an F-adapted 
process Q represent the (predictable) CSA value process of the contract, and an F-adapted 
process F stand for the value process of a CSA (cash) collateralization scheme. We denote 
by TT a real number meant to represent the wealth of the hedging portfolio of the bank in the 
financial interpretation. The close-out cash-flow R = R{tt) is in fact twofold, decomposing 
into a close-out cash-flow i?* from the bank to the counterparty, minus, in case of default of 
the bank, a cash-flow Rf = Rf (vr) from the funder to the bank (depending on vr). These two 
cash-flows are respectively derived from the algebraic debt x of bank to her counterparty 
and X(7r) of the bank to her funder, modeled at time r as 



X = Qr-Tr, X(7r) = -(vr - F^_). 



(2) 
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The close-out cash-flow, if r < T, is then modeled as Riir) = B} — lr=eR'^ (t^), where 

R^ = Tr + lr=e{pX^ - X-) - K=-e{PX- - X+) - ^e=-eX 
i?/(7r) = (l-r)X+(7r), 

in which p and /) stand for recovery rates between the two parties, and r stands for a recovery 
rate of the bank to her funder. The measurable exposure at default is defined in terms 
of R as 

^(Tr) := - Rin) = - + - p)x+ + (1 - r)X+(7r)) - (1 - p)lr=eX- , (4) 

where the second equality follows by an easy algebraic manipulation. 

We now consider the cash-flows required for funding the bank's position, meant in the 
sense of the contract and its hedging portfolio altogether. For simplicity we stick to the 
most common situation where the hedge is self- funded as swapped and/or traded via repo 
markets (see |Crl2a] ). The OIS rate rt is used as a reference for all other funding rates, 
which are thus defined in terms of corresponding bases to r^. Given such bases bt and bt 
related to the collateral posted and received by the bank, and Xt and related to external 
lending and borrowing, the funding coefficient 5t(vr) is defined by 

gtiTT) = {btTt - btTi) + Xt{7r- Tt)'' -Xtin- r^)' . (5) 

Then (rjvr -|- gt{'K))dt represents the bank's funding cost over {t,t + dt), depending on its 
wealth vr. 

Remark 2.3 A funding basis is typically interpreted as a combination of liquidity and/or 
credit risk, see [FT 11 a] and |CD12j . Collateral posted in foreign-currency and switching 
currency collateral optionalities can also be accounted for by suitable amendments to b and 
b, see [FSTinj and [PitT2] . 



2.3 BSDEs 



With the data ^ and g specified, the TVA process can be implicitly defined on [0, r] as the 
solution to the following BSDE, posed in integral form and over the random time interval 
[0,r]: For t G [0,f], 



/3flr<TCiPr - Qr-) + MPs - & s)ds 



(6) 



The reader is referred to jCrl2b| Proposition 2.1] for the derivation of the TVA BSDE Q. 
In this paper for simplicity of presentation we take ([g]) as the definition of the TVA. 

The practical conclusion of |Crl2b| is that one can even adopt a simpler "reduced, pre- 
default" perspective in which defaultability of the two parties only shows up through their 
default intensities, see Proposition 2.4 below and equation (3.8) in |Crl2bj . For t S [0,T] 
and TT G M let 



(7) 



Xt = Qt- Tt 

6(vr) = {Pt - Qt) +Pt{{l - Pt)xt + (1 - ^(vr - rt)-(7r)) -p,{l - p,)xT , 
in which 

p^ = p(r = 0|g,„), p^ = F{T = e\gr-). 

Note that in case of unilateral counterparty risk, we have 9 = oo and consequently, pr = 0, 
Pt = 1- 
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Proposition 2.4 (TVA Reduced- Form Representation) One has @ = @ on [0, r) 

and Qf = Itkt^,, where 



(3tet = Et / /3s{gs{Ps-es) + lsUPs-es))ds 



(8) 



forte[0,T]. 



Remark 2.5 This assumes that the data of ([s]) are in F, a mild condition which can always 
be met be passing to F-representatives (or pre-default values) of the original data. 

Remark 2.6 For r = 1 the exposure at default ^ does not depend on vr, and in case of a 
linear funding coefficient given as gt{P — t?) = QtiP) — ^^"(9, for some and A'^, the TVA 
equations ^ and ([s]) respectively boil down to the explicit representations 



(9) 

(10) 



Here the funding-adjusted discount factors are 

ft 



ft 



/3° = exp(- / (r, + A°)(is), /3° = exp(- / {7s + K)ds) 
Jo Jo 



and in ( 10 ) is given by 



= {Pt - Qt) + Pt{i - pt)xt - m - Pt)xt- 

On the numerical side such explicit representations allow one to estimate the corresponding 
"linear TVAs" by standard Monte Carlo loops provided that Pt and Qt can be computed 
explicitly. In general (for example as soon as r < 1), nonlinear TVA computations can 
only be done by more advanced schemes involving linearization [FTllbj . nonlinear regres- 
sion jCAC"'"10] or branching particles |HL12j . Deterministic schemes for the corresponding 
semilinear TVA PDEs can only be used in low dimension. 

In differential form the pre-default TVA BSDE ^ reads as follows (cf. |Crl2bL Defi- 
nition 3.1]) 

J Gt = 0, and for t £ [0, T] : 
\ -dQt = gtiPt - Qt)dt - djlt, 



(11) 



where fi is the F-martingale component of and 

gt{Pt -^)= gt{Pt -'&)+ itlm -^)- n'd. 
Remark 2.7 From ([T]), P satisfies the following F-BSDE 

J Pt = 0, and for t G [0, T] : 
I - dPt = dDt - nPtdt - dMt, 



(12) 



(13) 
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for some F-martingale M. Therefore, the following pre-default F-BSDE in 11 



follows from (11) and (13): 



Ut = 0, and for t G [0, T] : 

- dllt = dDt - {gt(^t) + nPt ) dt - duu 



p -e 



(14) 



where dut = dMt—d^f As the pre-default price BSDE ( 14 ) involves the contractual promised 



cash- flows dDt, it is less user- friendly than the pre-default TVA BSDE (11). This math 



ematical incentive comes on top of the financial justification recalled at the beginning of 
Section [2] for adopting a two-stage "clean price P minus TVA correction 0" approach to 
the counterparty risk and funding issues. 



2.3.1 Pre-default Markov Setup 

Assume 



gt{Pt-'&)=9{t.Xt,i 



(15) 



for some deterministic function g{t,x,9), and an M'^-valued F-Markov pre-default factor 
process X. Then 0j = Q{t,Xt), where the pre-default TVA pricing function Q{t,x) is the 
solution to the following pre-default pricing PDE: 



Q{T,x) =0, xeR'^ 

{dt + X) Q{t, x) + g{t, X, G(t, x)) = on [0, T) x R'^, 



(16) 



in which X stands for the infinitesimal generator of X. As mentioned in Remark 2.6, from 



the point of view of numerical solution, deterministic PDE schemes for (16) can be used 



provided the dimension of X is less than 3 or 4, otherwise simulation schemes for (11) are 
the only viable alternative. 



2.4 CVA, DVA, LVA and RC 

Plugging ([?]) into (12) and reordering terms yields 



gtiPt -{^)+rt^ = - jtPtil - MQt - Tt)' 
+ ltPt{{l-p){Qt-Tt)+ 
+ htVt -htV^ + \t{Pt-^-Tt) 
+ -it{Pt-'d- Qt) , 



\t{Pt-^-Tt 



(17) 



where the coefficient At := \t — 7tPt(l — t) of {Pt — "d — Tt)~ in the third line can be 
interpreted as an external borrowing basis net of credit spread. This coefficient represents 
the liquidity component of A. From the perspective of the bank, the four terms in this 
decomposition of the TVA can respectively be interpreted as a costly (non-algebraic, 
strict) credit value adjustment (CVA), a beneficial debt value adjustment (DVA), a liquidity 
funding benefit/cost (LVA), and a replacement benefit /cost (RC). In particular, the time-0 
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TVA can be represented as 



eo = -]E 



E 







l3Ml-p)Pt{Qt-Tt)-dt 



ptit{i- p)pt{Qt-rt)+dt 



E 



+ E 



(18) 



pt{btTt - ItT^ + XtiPt -Qt- rt)+ - XtiPt -Ot- Tty)dt 



Pat{Pt -Qt- Qt)dt 



The DVA and the ^tPt{^ — ^)iPt — Qt — Fj)" -component of the LVA can be considered as 
"deal facihtating" as they mcrease the TVA and therefore decrease the price (cost of the 
hedge) the bank can consider selhng the contract to her counterparty. Conversely, the CVA 
and the Xt{Pt — ©i — ^t)~ components of the LVA (for Xt positive) can be considered as 
"deal hindering" as they decrease the TVA and therefore increase the price (cost of the 
hedge) for the bank. The remaining terms can be interpreted likewise as "deal facilitating 
or hindering" depending on their sign, which is unspecified in general. 

3 CSA Specifications 



In the next subsections we detail various specifications of the general form (17) of g, de 



pending on the CSA data: the close-out valuation scheme Q, the collateralization scheme 
r and the collateral remuneration bases b and b. 

3.1 Clean CSA Recovery Scheme 



In case of a clean CSA recovery scheme Q = P, (17) rewrites as follows: 



UPt - 1?) + nt? = - -ftPtii - p){Pt - ^tr 
+ ^,p,((i-p)(p,-rt)+ 

+ btTf-btTj + \t{Pt-^-rt) 



(19) 



XtiPt -Ft 



Note rt on the left-hand side as opposed to in (17). In case of no collateralization, i.e. 
for r = 0, the right-hand-side of (19) reduces to 

-^tPtii - p)Pt' + 7tpt(i - p)Pt'' + [Pt - ^?)+ - ~Xt [Pt - ; (20) 

whereas in case of continuous collateralization with T = Q = P, it boils down to 

-btPt' + btPt+ + Xt^~ -Xt^+. (21) 

Remark 3.1 If A = A (case of equal external borrowing and lending liquidity bases), the 
TVA is linear for every collateralization scheme of the form Pt—Tt = et, for some exogenou^ 
residual exposure et, cf. Remark |2.6[ Setting 

9t = - iltPti^ -P) + At) + [jtPtil -P) + Xt) et + btTf - kV^, (22) 

'^Not depending on St, like with null or continuous collateralization. 



one ends up, similarly to (10), with 



(23) 



for the funding-adjusted discount factor 



'0 



(24) 



3.2 Pre-Default CSA Recovery Scheme 



In case of a pre-default CSA recovery scheme Q = H = P — @, (17) rewrites as follows 



+n^ = - [jtPtil -P) + At j {Pt-^- Tt) 

+ htPtii-p) + Xt){Pt-^-rty 
+ btTt-btrT. 



(25) 



In case of no coUateralization, i.e. for F = 0, the right-hand-side reduces to 



(7tPt(l -P) + At) {Pt - + iltPtil -P) + At) (Ft - ^)+; 



(26) 



whereas the continuous coUateralization with T = Q = P — Q yields 



bt{Pt-i})+ -btiPt-i^y 



(27) 



Remark 3.2 li b = b (case of equal collateral borrowing and lending liquidity bases), the 
TVA is linear for every coUateralization scheme of the form Pt — Qt — Tt = £t, for some 
exogenou^ residual exposure et, cf. Remark 2^ Setting 



9t=-{ m(l - P) + At ) Et + bt {Pt - St) + iltPtil -p) + At) et, (28) 



one ends up, again similarly to (10), with 



/3fet = ]Et f&,ds 



(29) 



for the funding-adjusted discount factor 



/3l' = exp(- / {r, + bs)ds). 
Jo 



(30) 



Remark 3.3 If 6 = 6 = 0, the two continuous coUateralization schemes of equations (21)) 



and (27) equally collapse to B = and H = P = T = Q. 



Not depending on Qt, like with continuous coUateralization. 
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3.3 Full Collateralization CSA 



Let us define the full collateralization CSA hy Q = T* where T* is given as the solution to 
the following F-BSDE: 



r^. = Pt, and for t G [0,r] : 

dr* - {{n + bt){r:)+ - {n + 6t)(r*)-) dt = dPt - nPtdt + d/l* 



(31) 



for some F-martingale /x*. Then Q := P — T* solves the pre-default TVA BSDE (11), or in 



other words 11 = F*. Indeed one has by (17), that for t G [0, T] 



g,(r*) + niPt-T* 



hence 



in + 6t)(r*)+ - in + hmr - nPt = gti^)- 



Therefore, the second line of (31) reads as 

-idp - dr;) = gti^;)dt - dji; = miPt - in - r*))dt - dji;, 

which, together with the fact that P — T* vanishes at T, means that P — T* satisfies the 



pre-default TVA BSDE Jllj). 

If 6 = 6, thenJtiUt) + nPt reduces to in + bt)Ilt in the BSDE ^ for the fully 
collateralized price 11 = F* . This BSDE is thus equivalent to the following explicit expression 
for B: 



IE. 



(32) 



where the funding-adjusted discount factor is defined by (30). In the special case b = 
6 = 0, we have T* = P which is the situation already considered in Remark |3.3[ This case, 
which yields G = and Yi = P = Q = T , justifies the status of formula ([T]) as the master 
clean valuation formula of a fully collateralized price at an OIS collateral funding rate rj. 
With such a fully collateralized CSA there is no need for pricing-and-hedging a (null) TVA. 
The problem boils down to the computation of a clean price P and a related hedge, see e.g. 
|Crl2b| for possible clean hedge specifications. 

3.4 Pure Funding 

In case 7 = 0, which is a no counterparty risk, pure funding issue case, the CSA value 



process Q plays no actual role. In particular the dt-coefficient of the BSDE (14) for 11 is 
given by 

gtiiit) + nPt = in + bt)Tt - in + bt)Ti + in + At)(n - Tt)+ - in + At)(n - r^)-, 

where A = A is a pure liquidity external borrowing basis. In case E = and A = A this 
results in the following explicit expression of H: 

^^Bi = IEj,^^^dZ)„ (33) 



with the funding-adjusted discount factor /3 defined in (24). In the special case A = A = 
one recovers the classical valuation formula ([T]) for 11 = P. 
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3.5 Asymmetrical TVA Approach 

In practice the bank can hardly hedge her jump-to-default and therefore cannot monetize 
and benefit from her default unless and before it actually happens). If one wants to ac- 
knowledge this, one can avoid to reckon any actual benefit of the bank from her own default 
by letting p = r = 1, cf. Q. Such an asymmetrical TVA approach, even though still 
bilateral, allows one to avoid many concerns of a general symmetrical TVA approach, such 
as the arbitrage issue that arises for r < 1, the hypothetical and paradoxical benefit of the 
bank at her own default time, and the puzzle for the bank of having to hedge her own 
jump-to-default risk in order to monetize this benefit before her default (see |Crl2aj and 



|Crl2bj ^'). Indeed in this case equation nlTn reduces to 



gtiPt -^) + rt^ = - - MQt - Tt)' 

+ btTf - btV^ + At (Pt - i9 - rt)+ - At (P - I? - Tt)- (34) 
+ jt{Pt-^- Qt) , 

where there is no beneficial debt valuation adjustment anymore (for p = 1 the second line of 



(17) vanishes), and where the borrowing funding basis is interpreted as a pure liquidity 
cost. ^ ^ 

Note that in case of a pre-default CSA recovery scheme Q = n_ = P — 0_, an 
asymmetrical (but still bilateral) TVA approach is equivalent to a unilateral TVA approach 
where the bank would simply disregard her own credit risk. One has in both cases 

gt{Pt -^) + = htVt - htV^ + \t{Pt-'&- rt)+ - {\t + 7iPt(l - 'P)) {Pt-'9- Tt)- , 

where jp is the intensity of default of the counterparty. Moreover, in an asymmetrical TVA 
approach with Q = n_, an inspection of the related equations in |Crl2bj shows that a 
perfect TVA hedge by the bank of an isolated default of her counterparty (obtained as the 
solution to the last equation in |Crl2bj ) in fact yields a perfect hedge of the TVA jump-to- 
default risk altogether (default of the counterparty and/or the bank). This holds at least 
provided that the hedging instrument which is used for that purpose (typically a clean CDS 
on the counterparty) does not jump in value at an isolated default time of the bank, a mild 
conditiorr] satisfied in most models. 



4 Clean Valuations 

In the numerical Section [5j we shall resort to two univariate short rate models, presented 



in Subsections 4.2 and |4.3[ for TVA computations on an interest rate swap. These are 



Markovian pre-default TVA models in the sense of Subsection 2.3.1, with factor process 
Xt = Tt. Our motivation for considering two different models is twofold. Firstly, we want to 
emphasize the fact that from an implementation point of view, the BSDE schemes that we 
use for TVA computations are quite model-independent (at least the backward nonlinear 
regression stage, after a forward simulation of the model in a first stage). Secondly, this 
allows one to assess the TVA model risk. 

Remark 4.1 The choice of interest rate derivatives for the illustrative purpose of this 
paper is not innocuous. The basic credit risk reduced-form methodology of this paper is 

In the notation of the concluding Subsection 4.4 of [Crl2bj . this is given 
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suitable for situations of reasonable dependence between the reference contract and the two 
parties (reasonable or unknown, e.g. the dependence between interest rates and credit is 
not liquidly priced on the market, see |BP08] ). For cases of strong dependence such as 
counterparty risk on credit derivatives, the additional tools of |Crl2cj are necessary. 

Talking about interest rate derivatives, one should also mention the systemic counter- 
party risk, referring to various significant spreads which emerged since August 2007 between 
quantities that were very similar before, like OIS swap rates and LIBOR swap rates of differ- 
ent tenors. Through its discounting implications, this systemic component of counterparty 
risk has impacted all derivative markets. This means that in the current market conditions, 
one should actually use multiple-curve clean value models of interest rate derivatives in the 
TVA computations (see for instance |CGN12) and the references therein). In order not to 
blur the main flow of argument we postpone this to a follow-up work. 



4.1 Products 

In the sequel we shall deal with the following interest rate derivatives: forward rate agree- 
ments (FRAs), IR swaps and caps, whose definitions we provide below. The latter are used 
in Section [5] for calibration purposes. The underlying rate for all these derivatives is the 
LIBOR rate. We work under the usual convention that the LIBOR rate is set in advance 



and the payments are made in arrears. As pointed out in Remark 4.1, we do not tackle 
here the multiple-curve issue. Thus, we use the classical definition of the forward LIBOR 
rate Lt{T, T + 6), fixed at time t <T, for the future time interval [T, T + 5]: 



6 \Bt{T + 6) 

where B.[T) denotes the time-t price of a zero coupon bond with maturity T. 



Definition 4.2 A forward rate agreement (FRA) is a financial contract which fixes the 
interest rate K which will be applied to a future time interval. Denote by T > the future 
inception date, by T + (5 the maturity of the contract, where 5 > 0, and by the notional 
amount. The payoff of the FRA at maturity T + 5 \s equal to 

pfra^j, j^^.r^^j, j^^^ ^ N5{Lt{T, T + 5)-K). 

The value at time t G [0, T] of the FRA is given by 

pf'-''{t;T,T + 5,K,N) = N{Bt{T) - KBtiT + S)), (35) 

where K = 1 + 6K. 



Definition 4.3 An interest rate (IR) swap is a financial contract between two parties to 
exchange one stream of future interest payments for another, based on a specified notional 
amount N. A fixed-for-fioating swap is a swap in which fixed payments are exchanged 
for floating payments linked to the LIBOR rate. Denote by Tq > the inception date, 
by Ti < • • • < T„, where Ti > Tq, a collection of the payment dates and by K the fixed 
rate. Under our sign convention recalled in Section [2] that the clean price values promised 
dividends dDt from the bank to her counterparty, the time-i clean price of the swap Pt for 
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the bank when it pays the floating rate (case of the so-called receiver swap for the bank) is 
given by 



Pt = P'^^it; n,Tn) = nI Bt{To) - Bt{Tn) - 5k-iBt{Tk) , 



(36) 



fc=i 



where t < Tq and S^-i = — T^-i- The swap rate Kf, i.e. the fixed rate K making the 
value of the swap at time t equal to zero is given by 



Kf 



BtjTo) - Bt(Tn) 
Y2=i h-iBt{Tk 



(37) 



The value of the swap from initiation onward, i.e. the time-t value, for Tq < t < T^, oi the 
swap is given by 



P'''\t;Ti,Tn) = N 



1 



B 



- K5k,-i Bt{Tk,) - Bt{Tn) -K J2 ^k-iBt{Tk) , 



fe=fet+i 



(38) 



where Tj^^ is the smallest (strictly) greater than t. If the bank pays the fixed rate in the 
swap (case of the so-called payer swap from the bank's perspective) , then the corresponding 
clean price Pt is given by Pt = -P^'"{t; Ti, r„). 

Definition 4.4 An interest rate cap (respectively floor) is a financial contract in which 

the buyer receives payments at the end of each period in which the interest rate exceeds 
(respectively falls below) a mutually agreed strike. The payment that the seller has to make 
covers exactly the difference between the strike K and the interest rate at the end of each 
period. Every cap (respectively floor) is a series of caplets (respectively floorlets). The 
payoff of a caplet with strike K and exercise date T, which is settled in arrears, is given by 

P^p\T-T, K) = 5{Lt{T,T + 5) - K)+. 

The time-i price of a caplet with strike K and maturity T is given by, with K = 1 + 6K, 



P'^\t; T, K) = 5Bt{T + 5)^^'^*' \{Lt{T, T + 6) - K)- 



£t 



= Bt{T + 5)^ 
= KBt{T)E 



pT+6 



Bt{T + S) 



- K 



K 



Bt{T + 5) 



= KM 



exp-i't'^'"^'^" 



K 



-Bt{T + 5) 



£t 



(39) 



The next-to-last equality is due to the fact that the payoff y ^^^^j^j^^-^ — Kj at time T + S 

kV = K - Bt{T + 5))^ at timer. The 



is equal to the payoff Bt {T + 6) . 
last equality is obtained by changing from the forward measure P"^ to the spot martingale 
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measure P, cf. [MROSj Definition 9.6.2]. The above equalities say that a caplet can be seen 
as a put option on a zero coupon bond. 

In the two models considered in the next subsections, the counterparty clean price P 



of an interest rate derivative satisfies, as required for (15), 

Pt = P{t,Xt) 

for all vanilla interest rate derivatives including IR swaps, caps/floors and swaptions. 



(40) 



4.2 Gaussian Vasicek short rate model 

In the Vasicek model the evolution of the short rate r is described by the following SDE 

drt = a{k - rt)dt + dWt , 

where a,k > and is a Brownian motion with volatility cr > on the filtered probability 
space (0,^r,F,P). The unique solution to this SDE is given by 



rt = roe""* + k{l - e"'^*) + T e-"(*-")(iW„^. 

Jo 



The zero coupon bond price Bt{T) in this model can be written as an exponential-affine 
function of the current level of the short rate r. One has 



where 



m™(t, T) = Q (l - e-"(^-*)) -T + ?j- 



(41) 
(42) 



with Roo = k — ^ , and 



n^a{t,T) := -e"^ j\-''''du = ^ [e-<^~^^ - l) . 

The clean price P for FRAs and interest rate swaps can be written as 

Pt = P{t,rt), tG[0,r], 



(43) 



inserting the expression (41 



for the bond price Bt{T) into equations (35), (36) and (38) 



from Definitions 14.21 and 14.3 

In particular, the time-t price, for Tq < t < T„, of the interest rate swap is given by 



{t,Tk)+n^a(t,Tk)rt 



(44) 



k=kt+l 



by (42) and (43). 



which follows from (38) and (41 ). In the above equation mya{t, Tfc) and nva{t, T^) are given 
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4.2.1 Caplet 



To price a caplet at time in the Vasicek model, one uses (39) with 

Brp{T + 6) = e'^vaiT,T+5)+n,a{{T,T+S)rT 



for m^a(T, r+ (5) and n„a(r, r + (5) given by (42) and (43). Combining this with Proposition 
11.3.1 and the formula on the bottom of page 354 in [ MR05j yields (recall K = 1 + 5K): 

P'P'iO; T,K) = Bo{TM-d^) - KBo{T + 5)^'(-d+), (45) 

where $ is the Gaussian distribution function and 

EVT 2 

with 

H^T := ^ (1 - ( 1 - e''^') ' . (46) 



2a3 

4.3 Levy Hull- White short rate model 

In this section we recall a one- dimensional Levy Hull- White model obtained within the HJM 
framework. Contrary to the Vasicek model, this model fits automatically the initial bond 
term structure Bq{T). 

As in Example 3.5 of |CGN12) . we consider the Levy Hull-White extended Vasicek 
model for the short rate r given by 

drt = a{K{t) - rt)dt + dZ^, (47) 

where a > and Z'' denotes a Levy process described below. Furthermore, 

<t) = fo{t) + -dtfoit) + A (- (e-"* (e-°* - 1)) ie-*, (48) 

where /o(i) = — 9tlog-Bo(t) and Tp^ denotes the cumulant function of Z'^] see |CGN12l 
Example 3.5] with the volatility specification (Ts{T) = e~"^'^~^\ < s < T, therein. 

In this paper we shall use an inverse Gaussian (IG) process = {Z^)t>o, which is a 
pure-jump, infinite activity, subordinator (nonnegative Levy process), providing an explicit 
control on the sign of the short rates (see |CGN12) ). The IG process is obtained from a 
standard Brownian motion W by setting 

Z^ = inf{s > ■.Ws + <is > t}, 

where ? > 0. Its Levy measure is given by 

1 ^2 



The distribution of is IG{-,t^). The cumulant function exists for all z G [— % 

2 

(actually for all z G (—00, ^] since is concentrated on (0, 00)) and is given by 



V',(z) = ?(l- Jl-2^). (49) 
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Similarly to the Gaussian Vasicek model, the bond price Bt{T) in the Levy Hull- White 
short rate model can be written as an exponential-affine function of the current level of the 
short rate r: 

Bt{T) = e'^Ut,T)+nUt,T)n ^ ^^g^ 

where 



mUt,T) := log(^^) -n(t,r) 



(51) 



a 



-a{T-s) 



a 



and 



niJt,T) := -e"* [ e-'^'^du = - 
Jt a 



-a(t-s) 



1 . 



ds 



(52) 



In the Levy Hull- White model the clean price P for FRAs and interest rate swaps can 
be written as 

Pt = P{t,rt), te[0,T], 



by combining the exponential-affine representation (50) of the bond price Bt{T) and Def- 
initions 4.2 and 4.3, In particular, the time-t price, for To < t < T^, of the swap is given 
by 



Pt 



N e 



k=kt+l 



(53) 



which follows from (38) and (50). In the above equations mie{t,Tk) and nie{t,Tk) are given 
by dsTl) and (pi). 



4.3.1 Caplet 



To calculate the price of a caplet at time in the Levy Hull- White model, one can replace 
B* with B, S* with S, A* with A, and insert S* = and A* = in Subsection 4.4 of 
|CGN12] . thus obtaining the time-0 price of the caplet 



P^P\{)-T,K) = Bq{T + 5)WF'^^' 
= i?o(T + 5)]EP"-'' 



1 



Bt{T + 5) 



K 



with 

Y := loe 



BojT) 
Bo{T + S) 



+ r {As{T + 6) - As{T))ds + r {T.,{T + 5)-^s{T))dZl 
Jo Jo 



where S^(t) = ^(l - e'^^*"*)) and As{t) = ^/'^(-^^(t)), for < s < The time-0 price of 
the caplet is now given by (cf. |CGN12| Proposition 4.5]) 



P^P'(0;r, K) 



2tt 



{iv - R){1 + iv - R) 



dv. 



(54) 



16 



for R > 1 such that My^^[R) < oo. The moment generating function My^^ of Y under 
the measure P'^"'"'' is provided by 



tT+5 



exp(^-^ i/j,{-J:s{T + 6))ds^ 
Bo{T) 



(55) 



X exp ( z \^og 

xexp( l' ^,((z- l)S,(T + 5) -zS,(r))ds 
for z € C such that the above expectation is finite. Alternatively, the time-0 price of the 



caplet can be computed as the following expectation (cf. formula (39)) 



exp-It^'rsds _ Bt{T + 6)^ 



(56) 



where r is given by (47) and Bt{T + 5) by (50). 



4.4 Numerics 

In Section [5] we shall present TVA computations on an interest rate swap with ten years 
maturity, where the bank exchanges the swap rate K against a floating LIBOR at the end 
of each year 1 to 10. In order to fairly assess the TVA model risk issue, this will be done 
in the Vasicek model and the Levy Hull- White model calibrated to the same data, in the 
sense that they share a common initial zero-bond term structure Bq (T) below, and produce 
the same price for the cap with payments at years 1 to 10 struck at K (hence there is the 
same level of Black implied volatility in both models at the strike level K) . Specifically, we 
set ro = 2% and the following Vasicek parameters: 

a = 0.25, k = 0.05, a = 0.004 

with related zero-coupon rates and discount factors denoted by Rq(T) and Bq(T) = exp{—T Rq(T)) . 



It follows from (41) and after some simple calculations. 



R*iT) 
fo(.T) 



Rr. 



) + 



a 



4a3r 



dT{TRl{T)) = k + e-'''^ {ro - k) 
_ae-^ (ro - A;) - - (1 - e-'^^) 



-aT 



An application of formula (37) at time yields for the corresponding swap rate the value 
K = 3.8859%. We choose a swap notional of = 310.136066$ so that the fixed leg of the 
swap is worth 100$ at inception. 

In the Levy Hull- White model we use a = a = 0.25 (same speed of mean-reversion 
as in the Vasicek model), an initial bond term-structure Bo[T) fitted to Bo{T) by using 
/o(T) = foiT) above in (48), and a value of <j = 17.570728 calibrated to the price in the 
above Vasicek model of the cap with payments at years 1 to 10 struck at K. The calibration 
is done by least square minimization based on the explicit formulas for caps in both models 
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reviewed in Subsections 4.2.1 and 4.3. 1[ After calibration the price of the cap in both models 
is 20.161$ (for the above notional yielding a value of 100$ for the fixed leg of the swap). 

The top panels of Figure [T] show 20 paths, expectations and 2.5/97.5-percentiles over 
10000 paths, simulated in the two models by an Euler scheme r for the short-rate r on a 
uniform time grid with 200 time steps over [0, lOjyrs. Note one does not see the jumps on 
the right panel because we used interpolation between the points so that one can identify 
better the twenty paths. 



O.OBSi , 1 , 1 1 1 1 , , 1 0.14 





Figure 1: 20 paths with 200 time points each of the short rate rt and of the clean price 
process Pt = P{t, rt) of the swap. Left: Vasicek model ; Right: LHW Model. 

The top panels of Figure [2] show the initial zero-coupon rates term structure Rq{T) and 
the corresponding forward curves /o(T), whilst the corresponding discount factors Bq{T) 
can be seen on the lower left panel. This displays an increasing term structure of interest 
rates, meaning that the bank will on average be out-of-the-money with a positive Pt = Pf^ 
in ( |38| ), or in-the-money with a negative Pt = —Pf^, depending on whether the bank pays 
floating (case of a receiver swap) or fixed rate (case of a payer swap) in the swap, see the 
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bottom panels of Figure [T] Note that the swap price processes have quite distinct profiles 
in the two models even though these are co-calibrated. The bottom right panel of Figure 
[2] shows the Levy Hull- White mean-reversion function K{t) corresponding to /q (T) through 





Figure 2: Top Left: Rl{T) ; Top Right: /o(r); Bottom Left: B^{T) ; Bottom Right: K{t). 



5 TVA Computations 

In this section we show with practical examples how the CVA, DVA, LVA and RC terms 
defined in Section [2] can be computed for various CSA specifications of the coefficient g in 
Section [3| The computation are done for an IR swap in the two models of Section [4j 

5.1 TVA Equations 



The generator X of the Gaussian Vasicek short rate r in the pre-default pricing PDE ( 16 ) 
is given by 

X@{t,r) =(a(A;-r))a,e + ^a^a^^e. 



Assuming (15), the corresponding TVA Markovian BSDE writes: G(r, ry) = 0, and for 
te[0,T] : 

-de{t, n) =g{t, ru e{t, rt))dt - d^t, rt))dW^ . (57) 
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Similarly, the generator X of the Levy Hull- White short rate r driven by an IG process 
is given by 

Xeit,r)={a{K{t)-r))drQ+ [ (G(t, r + e) - e(t, r))F,(de), (58) 

where stands for the Levy measure of Z^^^ Assuming (15), the corresponding TVA 
Markovian BSDE writes: e(T,rT) = 0, and for t G [0,r] : 

-de(t,rt) =?(t,n,e(t,n))(it- / (e(t,rt_+e)-e(t,rt_))iV^(dt,de), 

where A^"^ stands for the compensated jump measure of Z'' . 

5.2 BSDE Scheme 

Even though finding deterministic solutions of the corresponding PDEs would be possible 
in the above univariate setups, in this paper we nevertheless favor BSDE schemes, as they 
are more generic - in real-life higher-dimensional applications deterministic schemes cannot 



be used anymore. We s olve (57) and (58) by backward regression over the time-space grids 
generated in Subsection 4.4 see the top panels of Figure [ij We thus approximate 0t(ti;) in 



(57) and (58) by Q - on the corresponding time-space grid, where the time-index i runs from 



1 to n = 200 and the space-index j runs from 1 to m = 10^. Denoting by Gj = i^\)\<j<m 
the vector of TVA values on the space grid at time i, we have = 0, and then for every 
i = n — 1, ■ ■ ■ ,0 and j = 1, • • • ,m 

for the time-step h = ^ = 0.02y (one week, roughly). The conditional expectations in space 
at every time-step are computed by a g-nearest neighbor average non-parametric regression 
estimate (see, e.g., |HTF09] ). with = 5 in our numerical experiments below. 

5.3 Numerics 

We set the following TVA parameters: 7 = 10%, b = b = X = 1.5%, A = 4.5%, p = 
50%, p = 70% and we consider five possible CSA specifications in this order: 



(l^, p, p) = 


(40,40,40)%, 


Q 


= P, 


r 


= 


(1:, p, p) = 


(100,40,40)%, 


Q 


= P, 


r 


= 


(1^, p, p) = 


(100, 100,40)%, 


Q 


= P, 


r 


= 


(t, p, p) = 


(100,100,40)%, 


Q 


= n. 


r 


= 


(t^, P, p) = 


(100,40,40)%, 


Q 


= P, 


r 


= Q 



p. 

Note that under the first CSA specification, one has A = 4.5% — 0.6 x 0.5 x 10% = 1.5% = A, 
so this is a linear TVA special case of Remark |3.1[ where the TVA at time can be computed 
through a straight Monte Carlo simulation. 

Moreover, we shall study the TVA in the two co-calibrated Vasicek and Levy models, 
and for the receiver and payer swaps. We thus consider twenty cases (5 CSA specifications 
X 2 models x receiver versus payer swap). 



The integral converges in (581 under technical conditions stated in |CGN12j . 
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Table [T] shows the time-0 TVAs and the corresponding CVA/DVA/LVA/RC decom- 
positions (four terms on the right-hand side of (18)) in each of the twenty cases. For 



benchmarking the numerical BSDE results we display in Figure [3] the time-0 TVA BSDE 
value versus the TVA Monte Carlo mid- and 95%-lower and upper bounds in the first CSA 
specification (linear). In all four cases the BSDE time-0 value of the TVA is close to the 
middle of the confidence interval. 



1.36 
1.34 
1.32 
1.3 
1.28 



-1.87» 



-2.05 



■1.9. 



Figure 3: Time-0 linear TVA BSDE value versus Monte Carlo mid-value and 95%-bounds. 
From left to right: Receiver swap in the Vasicek model; Receiver swap in the Levy model; 
Payer swap in the Vasicek model; Payer swap in the Levy model. 

The numbers of Table [l] are fully consistent with the CVA/DVA/LVA/RC interpreta- 
tion of the four terms in the TVA decomposition on the right-hand side of (|18|), given the 



increasing term structure of the data discussed in Subsection |4.4[ For instance, a "high" 
DVA of 1.75 in the Vasicek model and 2.34 in the Levy model for the receiver swap in the 
first row of Table [T| is consistent with the fact that with an increasing term structure of 
rates, the bank is on average out-of-the-money on the receiver swap with a positive Pt = Pf^ 
(see Figure [T]). The CVA, on the contrary, is moderate, as it should be for a receiver swap 
in an increasing term structure of interest rates, and higher ( "more negative" ) in the Levy 
than in the Vasicek model (-0.90 versus -0.06). The numbers of Table [l] are not negligible 
at all in view of the initial value of 100$ of the fixed leg of the swap. In particular, the 
LVA terms are quite significant in case of the payer swap with r and/or p = 100%, see the 
corresponding terms in rows 3 and 4 in the two bottom parts of Table [T] The choice of r 
and p thus has tangible operational consequences, in regard of the "deal facilitating" ("deal 
hindering") interpretation of the positive (negative) TVA terms as explained at the end 



of Subsection 2.4, It is worthwhile noting that all this happens in a simplistic toy model 
of TVA, in which credit risk is independent from interest rates. These numbers could be 
even much higher (in absolute value) in a model accounting for potential wrong way risk 



dependence effects between interest rates and credit risk, see Remark 4.1 

Figures [4] and [5] (receiver swap in the Vasicek and Levy model, respectively) and Figures 
[g] and [T] (payer swap in the Vasicek and Levy model, respectively) show the "expected 
exposures" of the four right-hand side terms of the "local" TVA decompositions (17) with 
d replaced by 0^ therein. These exposures are computed as space- averages over 10^ paths 
as a function of time t. Each time-0 integrated term of the TVA in Table [T] corresponds 
to the surface under the corresponding curve in Figures |4] to [7] (with mappings between, 
respectively: Figure [4] and the upper left corner of Table [T| Figure [5] and the upper right 
corner of Table [l} Figure [6] and the lower left corner of Table [l| Figure [7] and the lower right 
corner of Table [1]). 

Finally, Figures [s] (receiver swap) and [9] (payer swap) show the TVA processes in the 
same format as the swap clean prices at the bottom of Figure [T] 
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Figure 4: Receiver swap in the Vasicek model. Columns: CVA/DVA/LVA/RC, the 4 wings 
of the TVA; Rows: 5 CSA Specifications. 
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Figure 5: Receiver swap in the Levy model. Columns: CVA/DVA/LVA/RC, the 4 wings of 
the TVA; Rows: 5 CSA Specifications. 
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Figure 6: Payer swap in the Vasicek model. Columns: CVA/DVA/LVA/RC, the 4 wings of 
the TVA; Rows: 5 CSA Specifications. 
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5 10 5 10 5 10 O 5 10 




5 10 5 10 5 10 O 5 10 




Figure 7: Payer swap in the Levy model. Columns: CVA/DVA/LVA/RC, the 4 wings of 
the TVA; Rows: 5 CSA Specifications. 




Figure 8: Receiver swap: 20 paths with 200 time points each of the TVA process Qt = 
Q{t,rt). Left: Vasicek model ; Right: LHW Model. Top to Bottom: 5 CSA Specifications. 
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Figure 9: Payer swap: 20 paths with 200 time points each of the TVA process @t = Q{t, 
Left: Vasicek model ; Right: LHW Model. Top to Bottom: 5 CSA Specifications. 
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CVA 


DVA 


T T r A 

LVA 


RC 


1A7 


-O.Ud 


1.75 


0.71 


-0.92 


1.34 


f\ C\f\ 

-0.90 


2.64: 


0.72 


-0.85 


1.40 


-0.06 


1.75 


0.64 


-0.91 


0.93 


-0.90 


2.34 


0.15 


-0.68 


0.40 


-0.06 


0.00 


0.76 


-0.29 


-0.45 


-0.90 


0.00 


0.32 


0.12 


O.OD 


r\ no 

-O.Oo 


n c\c\ 

0.00 


0.74 


r\ c\f\ 

0.00 


-0.43 


-0.76 


n nn 

o.uo 


0.32 


0.00 


0.43 


0.00 


0.00 


0.72 


-0.29 


0.44 


0.00 


0.00 


0.72 


-0.29 


TVA 


CVA 


DVA 


LVA 


RC 


TVA 


CVA 


DVA 


LVA 


RC 


-1.90 


-2.45 


0.04 


-0.68 


1.17 


-2.08 


-3.28 


0.64 


-0.66 


1.25 


-2.64 


-2.45 


0.04 


-1.92 


1.67 


-3.17 


-3.28 


0.64 


-2.41 


1.92 


-2.67 


-2.45 


0.00 


-1.92 


1.68 


-3.59 


-3.28 


0.00 


-2.38 


2.11 


-3.59 


-1.77 


0.00 


-1.83 


0.00 


-4.80 


-2.49 


0.00 


-2.26 


0.00 


-0.50 


0.00 


0.00 


-0.81 


0.31 


-0.51 


0.00 


0.00 


-0.81 


0.31 



Table 1: Time-0 TVA and its decomposition. Top: Receiver swap ; Bottom: Payer swap. 
Left: Vasicek model; Right: LHW Model. 

Conclusion 

In this paper, which is a numerical companion to |Crl2aj and [Crl2b| . we show on the 
standing example of an interest-rate swap how CVA, DVA, LVA and RC, the four "wings" 
(or pillars) of the TVA, can be computed for various CSA and model specifications. Positive 
terms such as the DVA (resp. negative terms such as the CVA) can be considered as 
"deal facilitating" (resp. "deal hindering") as they increase (resp. decrease) the TVA 
and therefore decrease (resp. increase) the price (cost of the hedge) for the bank. Beliefs 
regarding the tangibility of a benefit-at-own-default, which depends in reality on the ability 
to hedge and therefore monetize this benefit before the actual default, are controlled by the 
choice of the "own recovery-rate" parameters p and r. Larger p and r mean smaller DVA 
and LVA and therefore smaller TVA, which in principle means less deals (or a recognition 
of a higher cost of the deal). This is illustrated numerically in two alternative short rate 
models to emphasize the model-free feature of the numerical TVA computations through 
nonlinear regression BSDE schemes. The results show that the TVA model risk is under 
reasonable control for both co-calibrated models (models calibrated to the same initial term 
structure, but also with the same level of volatility as imposed through calibration to cap 
prices). We emphasize however that the latter observation applies to the "standard" case 
studied in this paper without dominant wrong-way and gap risks, two important features 
which will be dealt with in future work. 
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